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Abstract 

We use exact general solutions for the spatially flat FRW and 
the anisotropic Bianchi I cosmologies to show that generically un- 
coupled scalar fields cooperate to make inflation more probable, while 
the presence of several interacting fields hinders the occurrence of the 
phenomenon, in accordance with previous results based on particular 
power-law solutions. Similar conclusions are reached in the case of 
Bianchi VIo spacetimes, for power-law solutions which are proved to 
be attractors. 
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1 Introduction 



In many inflationary models the effective potential energy density of a scalar 
field is responsible for an epoch of accelerated inflationary expansion []J . Very 
often one assumes that inflation is driven by a scalar field of the Liouville 
form, i.e., a exponential potential, because this kind of potential arises in 
various higher-dimensional supergravity and superstring models |4|. ||, 
§0- 

Although there are many scalar fields in superstring theories, in the past 
it was often assumed that typically only one scalar field was responsible for 
the inflation, while those having higher exponents were quickly redshifted 
away. However, it has been found || that the so-called assisted inflation 
may occur when several scalar fields are present, even if each individual field 
is too steep to drive the inflation, provided that the fields are uncoupled and 
interact only through the geometry. On the other hand, if the fields interact 
directly with each other, the opposite effect may happen and the presence of 



cross couplings beteween fields may hinder inflation || [L(| |TT |. 

Assisted inflation has been mainly studied in power-law solutions (a oc t p ) 
for the spatially flat Friedmann- Robert son- Walker (FRW) cosmology, which 
can be shown M to be the late-time attractor for the evolution of this kind of 
model. (Recently, Green and Lidsey |12| discussed in the context of assisted 
inflation the late-time evolution in a general geometry.) 

The purpose of this work is to extend previous studies on multi-scalar 
field cosmologies in two directions: firstly, we will use general solutions (in- 
stead of the special power-law ones) to analyze if the presence of several 
fields generically helps or impedes inflation, and secondly, going beyond the 
aforementioned FRW cosmology, we will consider the anisotropic inhomoge- 
neous generalization given by Bianchi type I models. Thus in section EO we 



deal with n interacting scalar fields in a FRW spacetime and use the general 
solution to show that the larger is the number of interacting scalar fields the 
less likely is inflation. Section |2.2| makes plausible for more general scalar 
field potentials the results obtained in section |2.1| for exponential potentials. 
We do not know the solution for n non-interacting scalar fields in a FRW 
cosmology in the general case, but we use the discussion of the late-time at- 



tractor of section £]3] to restrict the analysis of uncoupled fields in section \[A 



to the particular case in which all fields are assumed to be equal, for which 
the general solution can be found. We show that non-interacting fields gener- 
ically cooperate to assist inflation. The density fluctuations corresponding 
to the last case are discussed in section |2.5| . General solutions of anisotropic 
Bianchi I cosmologies with interacting and uncoupled fields are used in the 
first part of section IO and subsection EO, respectively, to check that also in 
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those cases interacting fields make inflation more difficult while uncoupled 
fields assit it. The stability of power- law solutions is dicussed in [3.3| . Fi- 
nally in section |3] we turn to power-law solutions of the Bianchi VIo model, 
reaching again the same conclusions. 



2 The n-scalar field problem in a flat FRW 
spacetime 

In the following we will consider two kinds of problems in flat FRW space- 
times in which there are n homogeneous scalar fields driven by exponential 
potentials. First of all, we will assume that the scalar fields are interacting 
through a product of exponential potentials. Then we will consider the case 
in which the scalar fields are uncoupled because the potential is a sum of 
potentials involving a single field. 

2.1 The interacting n-scalar field problem in flat FRW 
spacetime 

The problem of n interacting homogeneous scalar fields, 0j, driven by a prod- 
uct of n exponential potentials Vi = Voie~ ki< ^ i minimally coupled to gravity 
in a flat Robertson- Walker spacetime, with metric 



ds 2 = -dt 2 + a 2 {t) dx 2 + dy 2 + dz 2 



is formulated by the system of equations 



3H 2 = 1<P 2 + V, (2) 



+ 3/70 - Vk = 0, (3) 
where H = a/a and the potential 

V{$) = Vae-**, (4) 

allows for interactions between the fields. V is the constant V01V02 • • • V 0n , 
and k = (ki, k 2 , k n ) is an n -component constant vector with respect to an 
orthonormal basis in the n -dimensional Euclidean internal space to which 
the vector = (0 X , 2 , 4> n ), built with the n minimally coupled scalar 
fields, also belongs. In the remaining of this section k, k 2 and cf 2 stand 

for k 



k ■ k, and • respectively. Potentials of this type are of interest 
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because they may be considered just as an approximation to a more complex 
potential. In fact, in higher-dimensional superstring theories, the scalar field 
is like one of the matter fields that contribute to the action and effective 
potential of the theory. Loop expansion [[|, or expansion in the number of 
interacting particles [ l3fl , of the action leads to a perturbative expression of 



the potential which is a summation of exponential terms || ||. From (||)-(^D, 
and discarding the trivial static metric solution H = 0, |0 2 + V = , we get 

H = -^ 2 . (5) 

Using this equation and the system (|2])-(|3]) one finds that the Klein-Gordon 
equations ([3]) have the first integrals 

$=Hk+^, (6) 

where c— (ci, C2, c n ) is an arbitrary vector integration constant. As (|6]) 
involves only geometrical quantities the Einstein-Klein-Gordon equations un- 
couple and the general parametric solution can be obtained [fHf] , JT5|, |T6| . 



One can easily verify that in this kind of models the Einstein-Klein- 
Gordon equations have power-law solutions 
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a oc tk 7 , (7) 

so that they inflate at all times when k 2 < 2. We will show below that the 
general solution has the same kind of behavior when the scale factor a is 
large enough. 

Inserting (||) in ([5]) we obtain the following second-order equation for the 
scale factor a(t) 

s + s m s + —^- s 2m+1 = 0, (tr/Tr/2), (8) 
4 cos"' a 

where 

— # — * 

COS (7 = — — , (9) 

ck 

m = —6/k 2 < 0, c = |c|, a and t have been replaced by the new variables s 
and r defined by 

a = s r = ckt cos cr, (10) 

and the dot stands for derivation with respect to r. Once s(r) is known one 
can compute, in principle, the scale factor a(r) from ([10] ), and the fields <p(t) 
from equations (^). 
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Equation (|8]) is a particular case of the second order nonlinear ordinary 
differential equation 

s + af(s)s + Pf(s) J f(s) ds + 1 f(s) = 0, (f f ) 

where f(s) is some real function and a, [3 and 7 are constant parameters. 
Depending on the values of k 2 our equation (H) corresponds to the next two 
cases: 



1. k 2 ^ 6 and 



f( 8 )=8 m , a = l, P = ^ L , 7 = 0. (12) 

4 cos a 



2. k 2 = 6 and 

f(s) = -, a = l, (3 = 0, 1= —L-. (13) 
s 4 cos^ cr 

Inserting the first integrals © in the Einstein equation ), we get a 
quadratic equation in the expansion rate H which has real solutions only 
when its discriminant is non-negative. This condition leads to 

P < 777^^71 < 7. ( 14 ) 



I _|_ 2a 6 V 



discarding the oscillatory solutions of equation flTTp, which would be obtained 
for a negative potential or /3 > 1/4. It follows that, if k 2 7^ 6, real solutions 
exist always for 

k 2 < k\ = -V. (15) 
sin a 

This sets a restriction on the integration constants and the exponent of the 
exponential-potential. 



The general solution of (11) can be obtained making the nonlocal trans- 



formation of variables [17] 



f(s)ds, v = jf(s)dr. (16) 

Under this transformation (M) becomes a linear inhomogeneous ordinary dif- 
ferential equation with constant coefficients 

z" + az +pz + j = 0, (17) 
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which for a > and (3 > is the equation of a damped harmonic oscillator 
in a constant external field. Here the ' indicates differentiation with respect 
to 7]. In our case the nonlocal transformation of variables fllEp is 



s m+l 



m + 1 



1 1 — i s di~ ■ 



for k 2 7^ 6, and 



Ins, rj = I —, (19) 



for /c 2 = 6. 

Taking a = 1 in ([17]) we get 



6ie A+,? + & 2 e A -", for A; 2 ^ 6; 

6i + &2 e_,? — 7^, for k 2 = 6; 



(20) 



where 



A ± = zi±£EM, (21) 



and &i, 62 ar e real integration constants. 
We see from (HOl) and (M) that 



Tit 

a (x z 3 (™+!) , (22) 

so that, while for k 2 > 6 the scale factor goes to as z — > 0, one has that 
for k 2 < 6 the scale factor goes to infinity when z — > 0. On account of 
equations (|20| ) one may assume that the latter occurs as 77 — > 770, f° r an 
appropriate value of r] . Let us now write r\ = r] + 6r] and expand (|20|) 
and ( P2| ) about z = and r] = r] , keeping only terms linear in 5i], so that 
5z = 5ri and 

a = 6z s(m+i) = 5r/ 3( m +i) . (23) 



Now, we see from (|10|) and (0) that 

<5?7 = a- 3 5r oc a _3 <5t. (24) 



If we use this in fl23|) , we get successively 

m rn 

a oc a™+i <5t 3 («+ 1 ) , (25) 

and 
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aoc5t"~=5t^. (26) 

We see then that if fc 2 < 2 the solution does inflate at least along some time 
interval. Now, since k 2 = k\ + k\ + ■ — h fc 2 , one concludes from our analysis 
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of the general solution that the larger is the number of the interacting scalar 
fields the less likely will be k 2 < 2 as well that inflation take place. On the 
other hand, for 6 < k 2 < k 2 , one has < (3 < 1/4 that implies that both 



A± < 0; thus Q18D and ( ^t|) tell us that the scale factor goes to infinity 
when z — > oo and rj — ► — oo (with no loss of generality we are taking both 
constants b%, hi > 0). By keeping only the dominant terms we have for very 
large negative rj that 

aoce"3(^+T). (27) 
This relation and (|4j), that is generally valid, yield 

— oc a~ 3 8t (28) 
a 

whose general solution is a oc t 1 ^ 3 . In this case the solution does not in- 
flate and it approaches to the free scalar field solution. These results are 
in accordance with those obtained with particular solutions by other au- 
thors H [nj [TT . 



It would be surprising if for k 2 = 6 the scale factor had a behavior 
drastically departing from that suggested by the previous analysis when 
k 2 E (Jfeg,6)U(6,2). In fact, we see from flT!]) and ([19]) that z = 3 log a 
so that, as a consequence of (p0|) , a goes to infinity as r] — > — 00. From (|20|) 
we see that log z ~ —rj, when rj — ► —00 and, by using 



8a 5z 6t 

oc — ~ —ori oc — , (29) 



a log a 
which can be written as 

— ~ a ~Ht log a. (30) 
a 

In the asymptotic regime in which a — >• 00 we have a 5a < a 2 8a/ log a < a 2 5a, 
so that from (|30"1) we get a 2 < t < a 3 and, finally, 

t 1/3 <a<t 1/2 , when a ^00. (31) 

We see that the solution does not inflate. 



2.2 More general potentials 

The results of the above section show that we can introduce an n -dimensional 
Euclidean internal vector space containing the n-component vector <fi = 
(01, 02, •••> <j>n) built with n minimally coupled scalar fields. Let us assume 
now that the potential has the general form 

V = V($), $ = (32) 
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In this case the Einstein equation (Q) remains unchanged, however the Klein- 
Gordon equation (|j) becomes 

<p + 3H<p + V'k = 0. (33) 

where the ' indicates derivatives with respect to the variable $. Since (|2]), (|32|) 
and ( PB| ) are invariant under rotations of the orthogonal axes in the Euclidean 
space above mentioned, we may choose the first axis of this internal space 
along the vector k. Then the Klein- Gordon equation splits into one equation 
for 0i = 

+ 3#0 + kV = 0, (34) 
and n — 1 free field Klein-Gordon equations for 2 = • • ■ = <p n = ip 

4> + 3Hip = 0. (35) 

From (|3~5D we obtain the first integral ip = Co/a 3 , where Cq is an arbitrary 
integration constant. Hence the original n-scalar field problem is equivalent 
to considering a self- interacting scalar field with stiff matter. In fact, the 
Einstein equation (0) now reads 

3H 2 = l<P 2 + V+^- 6 , (36) 

where is the sum of n — 1 positive defined integration constants. 

Taking into account that the scalar field depends only on t, its energy- 
momentum tensor may be written in the perfect fluid form 

T ik = (p<f, + p<t>)uiU k + p^gtk , (37) 

where 

= l^-V^). (38) 
The fluid interpretation of the scalar field has proven very useful in the study 



of the inflationary and Q- matter scenarios [I8|. In particular it leads to 
consider its equation of state p^ = (7^ — 1) p^. On the other hand, the state 
equation for fluid representing stiff matter is pf = (7/ — l)p/ with 7^ = 2. 
Because of the additivity of the stress-energy tensor it makes sense to consider 
an effective perfect fluid description with equation of state p = (7 — 1) p 
where p = Pf + p<t>, p = Pf + p<f> and 

7 = 7/P/ + 7 ^ (39) 
Pf + P4> 
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is the overall (i.e. effective) adiabatic index. For this effective perfect fluid 
the dynamical equations are 

a 2 , . 

3- = P 40 
a z 

and 

a 1 r . 
- = -g[p + 3p], (41) 

where p and p are the density and pressure of the effective perfect fluid. 
They involve the self-interacting scalar field and the n-l free scalar fields. 
Inflationary solutions occur when a > 0; this means that the expansion is 
dominated by a gravitationally repulsive stress that violates the strong energy 
condition p+3p < or equivalently 7 < 2/3 . When we impose this condition 
on (0), we obtain that 

V($) > 2 + % (42) 

Hence, when there is only one scalar field the solution inflates for V(0) > 
2 . However, the n interacting scalar fields produce a desassisted inflation 
because the solutions inflate only if ([42"D holds. Using the effective perfect 
fluid description we can generalize this analysis to an arbitrary potential 
V(<pi, 02> <f>n)- in this case we obtain that V(<px, 02, •••> <j>n) > 4>i + 01 + 

+ so that the interacting n scalar fields might make inflation more 

unlikely in a FRW spacetime. 



2.3 The n-scalar field attractor problem in the FRW 
spacetime 

Let us now assume that the n homogeneous scalar fields, 0j, in the FRW 
spacetime are driven by a general potential V = V(<j)i). In that case the 
Einstein-Klein-Gordon equations are 

1 n • 

3H 2 = - £ 0? + V, (43) 

i=l 

0. + 3 ij0. + V ^ = 0. (44) 
where stand for dV/dfa. From these equations we get 

1 n ■ 1 

#=-~£0;. (45) 
/ i=i 
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In order to investigate the stable scalar field configurations, we introduce 
the quantity 

u = jz—, (46) 

which reduces to u = 1 for the configuration <f>\ = 2 — • • • — 0n- Using 
d43f)-(]46D we find the differential equation for u: 

u = 2 nV ^-V. (47 ) 

n(f>l 

(In this section no summation convention applies to repeated Greek indexes.) 
If we further assume that the potential satisfies the condition 

V = nVjJ a , (48) 



w = 2^(lu - 1) (49) 



equation (W7f) becomes 



which has a fixed point solution: to = 1. Furthermore, the general solution 
of can be found using fl44]) : 

= 1 + (50) 

where c is an arbitrary integration constant. However, it is useful to express 
the last solution in terms of geometrical quantities with the aid of fl4"5p and 
). The final result is 



UJ 



nc \ 1 



1 + . (51) 

2a 6 H J V ; 



Evaluating (|51[) in the asymptotic regime, it can be easily shown that the 
particular solution u = 1 is an attractor for evolutions that behave asymp- 
totically as a oc t u with v > 1/3. 

This result strongly suggests that the special case in which all scalar 
fields are equal may be the late-time attractor of more general problems. In 
fact, this has been proved in || for potential (f|) with k\ — ■ • ■ — k n and 
Voi = • ■ ■ = Von, which guarantee assumption (fjg|). We did not use this fact 
in the previous section because we were able to use the general solution with 
no additional assumption on ki, Voi and 0j, but it will be useful to simplify 
somewhat the problem analyzed in the following section. 
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2.4 The non- interacting n-scalar field problem in the 
FRW spacetime 

Let us now assume that the n homogeneous scalar fields, fa, in the spacetime 
given by (0) do not interact directly, but are driven by a sum of n exponential 
potentials Vi = Poje~ In that case the Einstein-Klein- Gordon equations 
are 



3# 2 = E 



+ Vi 



+ 3Hfa - kiV t = 0. 



(52) 
(53) 



From now on we will consider a simplified problem in which k\ = &2 = 
• • • = k n = k and Vbi = V02 = • • • = Vo n = Vq. As discussed in the previous 
section, we can expect in this particular case that in the asymptotic evolution 
all scalar fields tend to a common limit. That this is actually the case has 
been proved in M. In consequence, in the remaining of this section we will 
take 0i = 02 = • • • = <Pn = 0, so that V\ = V% = ■ • • = V n = V = Voe~ fc * and 
equations (E2J)-(B3) reduce to 



3H 2 = -0 2 + nV, 
2 



+ 3H(p - kV = 0. 
One can easily get from (|54D-(|55D 



n 



H = --0 2 



and the first integral of the Klein-Gordon equation (|55| 

: k 



-H + 



n 



(54) 
(55) 



(56) 



(57) 



where c is an arbitrary integration constant. 

As in the previous section, the Einstein-Klein- Gordon equations have 
power-law solutions || 

aoct^, (58) 

so that this type of solutions inflates at all times when k 2 < 2n. We will 
now show that the general solution of the Einstein-Klein-Gordon equations 
( |54]) and (|55"D has the same kind of behavior when the scale factor a is large 
enough. 
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Inserting (0) in fl56|) we obtain the following second-order equation for 
the scale factor a(t) 

s + s m s + J s 2m+1 = 0, (59) 

where m = —6n/k 2 < 0, the dot means derivative with respect to r and we 
have used, instead of a and t, the new variables s and r defined by 

a = s s ; r = ckt. (60) 

We have again a particular case of equation ( |TT| ) with a = 1 and to 
find the general exact solution of (|59|) one has to consider the following two 
possibilities: 



1. k ^ 6n and 



Tfl A- 1 

/( S ) = s m , /? = — ^, 7 = 0. (61) 



2. fc 2 = 6n and 

/(*) = -, /3 = 0, 7 = 7- (62) 
s 4 

The general solution of ( |59"D can be obtained by performing the nonlocal 
transformation of variables given in (p^8|)-(p!9|) with the new value of m, which 
reduces (^9]) to the linear inhomogeneous ordinary differential equation with 
constant coefficients (|i~7l). 

If one now repeats the analysis of the previous section by systematically 
taking o = and using the new m, the same final results for z and a are 
obtained: just replace k 2 = k\ + • • • k 2 by k 2 /n in the quantities involved 
in ( |20"D and (p2|). One readily concludes that, when a — > oo, these models 
inflate if k 2 < 2n, so that the fields cooperate to make inflation more likely 
in the so-called "assisted inflation" , which was first discussed — but only for 
power-law solutions — in ||. 

This result could have been anticipated since the present case is included 
in the mathematical problem set in section |2J] by equations (|8|) and @ by 
just taking cos cr = ±1 (that corresponds to the simplified problem under 
consideration) and m = —6n/k 2 . 



2.5 Density fluctuations 

The fact that the contributions of the density fluctuations differ significantly 
in different inflationary universe models, has motivated a detailed study of 
all the alternatives. In this context, it is interesting to derive the spectral 
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indices for the perturbations that would be created during the periods of 
inflation described by the solutions given in the last section. It is well known 
that, for multi-scalar field models, the spectrum of the curvature perturbation 
reads 



19 



Ps= — 



dN dN 



2ix) 



(63) 



where N is the number of e-foldings of inflationary expansion remaining, 
and there is a summation over i and j. In the case considered in the previous 
section, where all the scalar fields are equal, (|63"D yields M 



„ , r , (H\ 2 1 H 2 



(64) 



J aH=k 



where H and <f> have to be evaluated at the time when the wave number of 
interest k leaves the horizon during inflation. Also in this case, the spectral 
index ns(k) defined as 

d lnP 5 



n s (k) = 1 + 



d Ink 



is given by JT9 



ns 



"4- 



(65) 



The availability of exact solutions allows us to express the relevant quanti- 
ties as functions of the variable r\ introduced in equation (H^). The scale 
factor reads (see fl22j) , which is satisfied for uncoupled scalar fields with 
m = —6n/k 2 , as one can easily see) 



a(rj) = [(m + 1) z(r])} 3 



(66) 



where z(rj) is given by (pD|)— (2T]). The spectrum Ps(k) obtained from ( 64]) 
is shown in Figure 1, for k = 3, c = —1, b\ = —1, 62 = 0.001 and different 
values of n. Here, inflation is possible for n > 5, and one can see that the 
peak of the spectral distribution moves towards the high frequency region as 
n increases. The corresponding spectral index ns is shown in Figure 2. From 
(j65|) and the general solutions (|20"D, (|21|) and (|66D, it can be shown that the 
value of ns —>■ 1 in the asymptotic region a > 1 as the number n of present 
fields increases. This feature is exemplified in Figure 2, where we can see how 
the larger is the value of n, the closer is the spectrum to the scale invariance 
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3 The n-scalar field problem in Bianchi type 
I models 



Now we turn to the general Bianchi type I model with n homogeneous scalar 
fields driven by exponential potentials. As we proceeded in the case of FRW 
spacetimes, we will first assume that the scalar fields are interacting through 
a product of exponential potentials, an then we will consider the case in which 
the scalar fields are uncoupled because the potential is a sum of potentials 
involving a single field. 

3.1 The interacting n-scalar field problem in the aniso- 
tropic Bianchi type I model 

The general Bianchi type I model is the anisotropic generalization of the 
spatially flat FRW universe expanding differently in the x, y, and z directions. 
In the usual synchronous form its line element is given by 

ds 2 = -dT 2 + a\(T) dx 2 + a\{T) dy 2 + a 2 3 (T) dz 2 . (67) 

For convenience we use the semiconformal coordinates 
dT 

dt = — , er = o|, G = a%a2, e p = — , (68) 
a 3 " a 2 



to cast the metric (|67j) into the form 

^ = e f(t) (_ dt 2 + ^ + ^(t) dx l + e -p(t) dy 2^ (6Q) 

We first consider, as in section |2.1| , n scalar fields <pi interacting directly 
through the exponential potential (|4j). The problem of n interacting homo- 
geneous scalar fields, 0«, driven by a product of n exponential potentials 
Vi = Voje~ fci ^% minimally coupled to gravity in the Bianchi I spacetime (|69D, 
is formulated by the following system of Einstein-Klein-Gordon equations 

a 

P=g. (70) 

i4(i) 2 -i^? 2 -^ < 72 > 

$+^-e f Vk = } (73) 
G 
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where a is an arbitrary integration constant. It can be easily seen that the 
vector 

Gk m , 

* = o2 + g' (74) 

(where m is a n-dimensional vector whose components are integration con- 
stants) is a first integral of the Klein-Gordon equation set (|73|). Inserting 
(|70D in equation fl74D the general solution of the Klein-Gordon equations is 
found: 

f=$o+p— + -\nG, (75) 
a 2 

where O is an arbitrary constant vector. Equations (|7C] )-([72|) along with 
(|75| ) uncouple and their solutions can be obtained if one is able to solve the 
following third-order equation for G 

GG 2 - GGG + Q " GG 2 + (m 2 + j\ G = 0. (76) 

Once G(t) is known, in principle one can compute p(t) and <p(t) from equa- 
tions (|70|) and ([75]) , respectively; /(t) is then obtained from (|7I]) . 

The Einstein-Klein-Gordon equations admit power-law solutions, G = t a , 
but they happen to be isotropic and, thus, equal to the ones discussed in 
section ^7T| . Thus we will analyze, instead, the general solution of (|76|) . 

Equation ([Tf]) has the first integral 

G^ + (K _i)G+^L- = c, (77) 
G G 

where G is an arbitrary constant and 

P 1 n 2 
K=---, M 2 = m 2 + -. (78) 
4 2' 2 v ; 

If instead of t and G we use the new variables z and r defined, for G 7^ 0, 

in 

G = z l l K , t=-l (79) 



^ + + ^^- 2 /^ = 0, (80) 



then, equation (|77|) becomes 

KM 2 
~~C 2 

where a prime denotes the derivative with respect to r. This equation is, once 
more, a particular case of flTTD and can be linearized by using the non-local 
transformation (|i~6|) , which in this case is 

y = / Z -VK dz = K z } rj = / Z -VK dr = --p (81) 



K-l' ' J a 
15 



for K ^ 1, and 



y = J z 1 dz = In z, 77 = J ' z 1 dr = — —p, (82 



for K — 1. If we take 

C 2 

for K ^ 1 and 



M 2 

/3 = (i^-l) 7 = , (83) 



M 2 

/3 = 0, 7 = £2-, (84) 

for .K" = 1, equation fl80|) reduces to two particular cases of equation ([H]) for 
a — 1. The trivial solution of this equation gives the implicit general solution 
of (76) which can be written, for arbitrary a, M and non-vanishing C, as 



G 

for K ^ 1, and as 



(de Ar? + C 2 e~ Xv )] (85) 



G = C x e-^ c ^ " (86) 

for K — 1. C\ and C2 are integration constants and A = \/l — 4/3/2. 

To check whether a model inflates, we will look at the sign of the decel- 
eration parameter q = —9~ 2 [39 + 9 2 ), where 9 = u a a is the expansion and 



9 = 9^ a u a , u a being the four-velocity of the cosmic fluid. Since in this case 
we are dealing with comoving coordinates, u a = (e - -^ 2 , 0, 0, 0^, one can see 
that, apart from a positive factor, the deceleration is 

q oc 9KG 2 - 3(2km + C) + (km - Cf + 9M. (87) 

We see from ( j35|) that when K < 1 (i.e., when k 2 < 6) G blows up 
for some value 77 = r] = ^ log (— C 2 /Ci), provided that CiC 2 < 0. If we 
expand (|85"D around this value (r/ = r] + 5rj) we get 

Got6r)T^, (88) 



and from @ and (|8T|) 



so that 



dr C , , 
V = ~G = ~G (89) 

d = ~^ G ' (X K^l 67] ' 1 - (90) 
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In consequence, when G —* oo the deceleration parameter flSTj ) is 

ff< x9Jf^^i, (when 5^0), (91) 

and there is inflation if K < 0, i.e., if /c 2 = k\ + k\ + - ■ - + k 2 < 2. We conclude 
that in these anisotropic universes also a greater number of interacting scalar 
fields makes inflation less likely. 



3.2 The non-interacting n-scalar field problem in the 
Bianchi type I model 



We will now assume that the n homogeneous scalar fields 0, in the metric ([67]) 
do not interact directly, but are driven by a sum of n exponential potentials 
Vi = Voje~ fe ^. To simplify the task of finding exact solutions of the Einstein- 
Klein-Gordon equations, we will further assume that k\ — &2 = • • • = k n = k, 
(f) 1 = (f) 2 = ... = (f) n = (f) a ndV 1 = V2 = --- = V n = V = V e~ k<t> , so that 
aforementioned equations can be written as 

a 

P=g, (92) 

<j> + - Fi <j>-ke f V = a, (95) 
G 

where a is an arbitrary integration constant. It is easy to check that 

■ k G m , . 

^TnG + G (96) 
is a first integral of the Klein-Gordon equation ( ]95|) , in terms of the new 
integration constant m. Inserting fl92f) in equation ( ^6|) the general solution 
of the Klein-Gordon equations is found: 

(j) — p — In Cr, (97) 

a 2n 



where 0o is an arbitrary constant. The equations (92)— ( 0) along with 
(|97|) uncouple and their solutions can be obtained if one is able to solve the 
following third-order equation for G 

GG 2 - GGG +(--!pj GG 2 +( m 2 + ^pjG = 0. (98) 
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Since this equation is the same as ( |76| ) once one replaces k 2 = k\ + k\ + 
■ • • + k 2 by k 2 /n, we may repeat the calculations of the previous section 
to reach the opposite conclusion: if several non-interacting scalar fields are 
present, they will cooperate to "assist" the inflation, which will be more likely 
and occurs for k 2 < 2n. 



3.3 Stability of power-law solutions in Bianchi type I 
model 

For many purposes it is interesting to investigate the stability of the solutions 
of (0). In particular, we hope that the solution representing an accelerated 
expansion of the universe, and the solutions that correspond to the assisted 
inflation, be stable. To this end we introduce the variable 



n 



h 



where h — ^, in equation (|76|) : 



(99) 



+ 



VL + K 



M 2 
hHT 2 



(n + i)h = o. 



(100) 



This equation has the fixed point solution Q = — 1. Note that equation ( |100| ) 
has also the fixed point solution Q = —K if G — > oo asymptotically. The 
corresponding asymptotic limits of these solutions can be obtained by solving 
( |9"P| ) for them. The final result is G oc t and G oc t x ' K respectively. Let us 
investigate the stability of these solutions when G blows up. 



From (fog) it is 

easy to see that Q = — 1 is unstable because expanding the solutions about 
it, Q = — 1 + e with e < 1, the sign of e depend on the slope of potential 
and the initial conditions. In fact, the corresponding solution G oc t does not 
satisfies Einstein equations (cf. (|93|)) and was introduced when multiplying 
(|94l) with G 2 G to obtain (p8|). On the other hand, the asymptotic solution 
Q = —K is stable because the dynamical equation for the perturbation e 



K 



Kt 



:ior 



near the attractor indicates that e decreases for K < 1. In particular, the 
inflationary solutions, that occur for K < 0, are stable. 
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4 The n-scalar field problem in a Bianchi VIo 
model 



The Bianchi VIo model can be written as follows: 

ds 2 = e m (-dt 2 + dz 2 ) + G(t) (e z dx 2 + e~ z dy 2 ) . (102) 

We first consider, as in previous sections, n scalar fields <pi interacting 
through the exponential potential (|j). The corresponding Einstein-Klein- 
Gordon equations are 

e' = (103) 

i-Ki) 2 "i /+ ^-^ (104) 

G ' 

<p + -(p-e f Vk = 0. (105) 
As formerly, one can check that the vector 

*=G2 + G' (106) 

where m is a n-dimensional arbitrary constant vector, is a first integral of 
the Klein- Gordon equation set (|105|) . By using this result and the value of 
/ one obtains from (|103|), we get 



(\ k 2 \ 1 
GG 2 - GGG + I - - — GG 2 + -GG 2 + m 2 G = 0. (107) 

The solutions of this equation are not known. However, investigating the 
stability of its fixed points, the asymptotic behavior of the general solution 
can be obtained in a simple way. In order to see whether assisted inflation 
works in Bianchi type VIo metrics it is sufficient to analyze the special case 
m 2 = 0. In terms of the variable Q defined by (P§D, equation ( |1 U7| ) becomes 



(Q + l)h = 0. (10* 



This equation has three fixed points: Qi = —1, ^2 = — K if h — > oo 
asymptotically, and fi 3 = 0, which correspond to G oc t, G oc t x l K and 
G oc e t//v/ ^ respectively. Now, we investigate the stability of these solutions 
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when G blows up. Expanding the solution about fixed points, that is, making 
Q = ^1,2,3 + e with e < 1, we get 



i = ^« (109) 

for Qi, which shows it is unstable, as in the case discussed in the previ- 
ous section. On the other hand, one obtains equation ( |101| ) for the linear 
approximation around Q2, and 




(110) 

in the case of f^. We conclude that the asymptotic solution Q2 is stable for 
K < 0, which means k 2 < 2 (for this set of potential slopes we have inflation), 
and is stable for K > 0. Note that G oc t x l K is only an asymptotic solution 



of ( |107| ), with m = 0; however, it acts as an attractor for all solutions that 
are close to it. 

In this special case in which m 2 = 0, one can readily see that the the 
deceleration parameter for the solutions G oc t x l K is (after recovering the 
implicit absolute value around t) 

^i^r 2 - 1 '*, cm) 

which is negative for — 1/2 < K < 0. Again, the more interacting scalar fields 
the less likely is inflation, which ensues when k 2 = k\ + k\ + • ■ ■ + k 2 < 2. 
If one assumes now that the scalar fields do not interact directly and one 



sets all the fields equal, as in sections |2]J and |372| , it is easily seen that the 
attracting power-law solutions inflate when k 2 < 2n, so that the scalar fields 
cooperate to assist inflation. 



5 Conclusions 

We have studied the effects of the appearance of more than one scalar field 
both in FRW spacetimes and in anisotropic Bianchi I cosmologies. Instead 
of using the important but particular power-law solutions, we have taken 
advantage of the general solution to analyze the generic behavior in FRW. 
In Bianchi I the power-law solutions are isotropic and, thus, very particular, 
so that the use of the general solution is even more illustrative. In all cases 
we have found, in agreement with calculations made by other authors with 
power-law solutions in FRW, that the existence of more than one scalar 
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field assists inflation provided that they are uncoupled and interact only 
through expansion. Also, in this case, the spectrum of density perturbations 
becomes closer to scale invariance as the number of fields increases. If, on the 
contrary, the fields interact directly with each other, inflation is less likely 
to occur. The same behavior has been obtained in Bianchi VIo universes, 
but (not having available the general solution) only for power-law solutions, 
which have been shown to be attractors when inflation arises. These results 
reinforce our belief that the presence of several uncoupled scalar fields in more 
general cosmologies fosters inflation, but that, in contradistinction, mutually 
interacting scalar fields tend to hinder the inflationary process. 
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Figures 




Figure 1. Spectrum of the curvature perturbations for k — 3, c = — 1, 
bi — — 1 and 6 2 = 0.001. As n increases (inflation is possible when n > 5) 
the spectral distribution peak is shifted to high frequencies. 
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Figure 2. The spectral index n${k) for the general solution approaches 1 in 
the asymptotic region a>las the number of fields increases. 
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